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Abstract. - We study the field dependence of the quasi-particle density of states, the thermo- 
dynamics and the transport properties in the vortex state of d-wave superconductors when a 
magnetic field is applied perpendicular to the conducting plane, specially for the low field and the 
low temperature compared to the upper critical field and transition temperature, respectively, 
H/Hc2 <C 1 and T/T c <C f. Both the superfluid density and the spin susceptibility exhibit 
the characteristic \/if-field dependence, while the nuclear spin lattice relaxation rate T^ 1 and 
the thermal conductivity are linear in field H. With increasing temperature, these quantities 
exhibit the scaling behavior in T/yH. The present theory applies to 2D /-wave superconductor 
as well; a possible candidate of the superconductivity in Sr2Ru04. 



Introduction. - Single crystals of high-T c cuprate superconductors like YBCO, Bi2212, 
etc. appear to provide the most useful testing ground for properties of unconventional 
superconductors m. In particular d-wave superconductivity has been established both in the 
hole-doped and the electron-doped high-T c cuprates ||, |I| . Due to the nodal structure in 
the d-wave order parameter, the specific heatH, the spin susceptibility and the superfluid 
density H have been predicted to behave like \/H in the vortex state where H is the magnetic 
field. Indeed the \f~H dependence of the specific heat in the vortex state of YBCO has been 
established |7|, pi |j| . On the other hand, the \f~H dependence of the superfluid density appears 
to have not been seen in spite of an elaborate muon spin rotation experiment [1C]. 

In the meantime, the calculational technique greatly improved [O, [l^, [l3| 14]. Therefore 
we can study the effect of a magnetic field within the semi-classical approximation almost 
analytically within the weak-coupling theory of d-wave superconductivity. Most of earlier 
works considered only the spherical Fermi surface [O, |l2|, [l3| . Here we assume that the Fermi 
surface of these d-wave superconductors is a weakly modulated cylinder [H. Then the present 
model should be also applicable to the recently discovered d-wave superconductors in k-(ET)2 
salts |l6|, 17 . The present model applies as well to the superconductivity in Sr2Ru04, if 
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the superconductivity is one of 2D /-wave states (i.e. A(k) ~ e ±I ^ cos(2</>), e ±l< ^ sin(2(/>) or 
e ±# cos(cfe3), where <f> is the angle between the quasi-particle momentum and a-axis in the 
plane and the quasi-particle momentum in the c-direction and c is the distance between 
layers of conducting plane in the c-direction.) which are considered by Hasegawa et al. [p"8f , 
The specific heat |l9| . NMR |2(J and the magnetic penetration depth measurement show 
clearly the presence of the nodal structure p2[ . The thermodynamics and the planar transport 
of these /-wave states are exactly same as the ones in d-wave superconductors when the field 
is applied perpendicular to the conducting plane. 

The object of the present paper is to study the quasi-particle density of states in the 
vortex state of d-wave superconductors when a magnetic field is applied perpendicular to the 
conducting plane, specially for the low field and the low temperature compared to the upper 
critical field and transition temperature, respectively, H/H C 2 <C 1 and T/T c <C 1. Then 
making use of the density of states, we calculate the specific heat, the spin susceptibility, the 
superfluid density and the nuclear spin lattice relaxation rate T^ -1 in NMR. Also for clarity, we 
limit ourselves to the superclean limit where the effect of impurity scattering is negligible. For 
E/A -C 1 the quasi-particle density of states is a simple function of E/e ~ E/y/H where E is 
the energy of a quasi-particle, A is the d-wave superconducting order parameter, e = vy/eH/2, 
e electron charge and v the Fermi velocity in the conducting plane. Then in the limit of T — > 0, 
the specific heat, the spin susceptibility and the change in the superfluid density behave like 
VH, while Tj" 1 behaves like H. Indeed, the y/~H dependence of the spin susceptibility and the 
H linear dependence of T^ 1 , are recently observed by NMR in the vortex state of underdoped 
YBCO |2^| . In addition, the thermal conductivity both in the superclean limit and the clean 
limit will be briefly discussed. 

For T ^ K, all these quantities exhibit scaling behavior as first discussed by Simon and 
Lee |24|. Actually the present model gives the explicit scaling functions, which should be 
readily accessible experimentally. 

Quasi-particle density of states. - Following the semi-classical approximation by Volovik[^| 
the quasi-particle density of states in the the vortex state of d-wave superconductor is given 

by 

N(E,H)/N ~±(\E\V\vq\) (I) 

where |-E|V|v-q| means the bigger one among \E\ and |v-q|. Also we assumed \E\, |v-q| <C A. 

Here |v ■ q| is the Doppler shift pa associated the pair momentum 2q and the Fermi velocity 

v 

v. In a magnetic field H II c, the Doppler shift is given as v • q = — cos</> where r is the 

2r 

distance from the center of the vortex and <p is the angle between v and q. Although v is 
parallel to one of nodal lines at low temperature, <fi runs from to 2n. Finally (. . .) of the 
Eq.(l) means the spatial average over r and <f>. This average is carried out over a unit cell of 
a square vortex lattice characteristic to <i-wave superconductors || d la Wigner Seitz. 



rdr ... (2) 



where 2d is the distance between vortices, d — l/Ve~H Q. Then Eq.(l) reduces to 

»,-/th t t\ / t\ t 4 e .E. 2v riT VH ,2 E . 

N(E,H)/N = --rg(-) s — J—l—g(-J-°—) (3) 



7T A' 



e 



with the scaling function g(s), 
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Fig. 1. - The scaling function g(s) in Eq.(4) is shown as a function of s = E/e. 

Fig. 2. - The specific heat data (O) by Nishizaki et al. [19] is fitted with the 
vH law. 



l(s) -i W + for S =f >1 

" ,_ l 171^ ^ +^(1 + 2s 2 ) sin" 1 s fors = £/e<l W 



and e = -VeH = — V~H and <& is a quantum of flux (~ 2.07 x 10 -11 T-cm 2 ). For E = 
2 2 V $ v ; 

0, we obtain N(0,H)/N = since g(s) = 1 + is 2 for s < 1. Also g(s) ~ f s for s > 1 

which means N(E, H)/N$ ~ i?/A for the high energy excitations E e. The scaling function 

E /_ 2 /*iT E 



g(s) versus s = y( = w ^7Jf) ^ s s h° wrL i n Fig- 1- A function similar to Eq.(4) has been 
obtained for a spherical Fermi surface in [Q . 

Thermodynamics. - Making use of the density of states, i.e. Eq.(3), we can determine the 
thermodynamic quantities at the low temperatures, T -C T c . First, the specific heat C S (T, H) 
is given by 

C<T,H )h „T - ^^Msec^) 



where 



3 f Z" 00 p 
f(T/e) = ^(-) 3 jf d ss 2 5 ( s )sech 2 ( — s) (6) 



and 7„T is the specific heat in the normal state. We have f(T/e) — > 1 for T/e < 1 and 
27<f (3) T 

/(T/e) = — r^- — for T/e > 1. With this limiting behavior of /(T/e), the specific heat takes 
An e 

C S /7„T = ^ 4- = — at low temperature T < e and C s /7 n T = for the 

7T ZA 7T W ZA 7T ^\ 



H = 0, respectively. In the limit of T — ► 0, the v i? dependence of the specific heat in YBCO 
has been discussed in @, H When we parameterized C s /T — A c \f~H and the specific heat 

in the absence of the field C 8 (T,Q)/T = aT, we can deduce v = 27 ^ J^-— ~ 2.28 x 10 6 

27r y $0 cv 

cm/s and 1.76 x 10 6 cm/s fro m M and ||, respectively. We may introduce an adjustable 
parameter a in front of v as in ||12|| . Then if we assume the Fermi velocity v ~ 10 7 cm/s, a is 
about 0.2 which is the same as deduced by Chiao et al. |2G|] . 

In the case of Sr 2 Ru04 the recent specific heat data by Nishizaki et a/.|ljj exhibits clearly 
the \H dependence as shown in Fig. 2. The deviation from the yH law for H < 0.01T is 
most likely due to disorders. 

Also Wang et al. [^7j studied the scaling behavior of the specific heat of YBCO 

C(T,H)-C(T,0) _ T T 27C(3) T 

Ae = ^7 ) = /( 7 )_ ^r7 
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Fig. 3. - The scaling function F(T/e) of specific heat in Eq.(7), the scaling function 7(T/e) of the 
spin susceptibility (or the superfluid density) in Eq.(lO) and the scaling function G(T/e) of the nuclear 
spin lattice relaxation rate in Eq.(13) are shown altogether as a function of T/e. 



|. ( l n2 )(^ + ... for|»l 



We show in Fig.3 T(^) versus — . The present model appears to describe the scaling behavior 
in YBCO |27j very well. Clearly a similar study of the scaling relation in the vortex state of 
Sr 2 Ru04 will be of great interest. 

Similarly, both the spin susceptibility x an d the superfluid density p s should exhibit the 
scaling behavior. 

= -j-HT/,) (8) 

TT A 

and 

h(T/e) = -- J d S5 ( s )sech 2 ( — s) (9) 
Here \n is the spin susceptibility in the normal state. We have h(T/e) — > 1 for T/e <C 1 and 

7T T 

h(T/e) = — \n2— for T/e > 1. With this limiting behavior of h(T/e) both the susceptibility 
and superfluid density take ~ 2s i n th e limit of T — > 0, and ~ 2 (In 2)-^ in the limit of 



T > e As already mentioned in the introduction the v if dependence of the susceptibility has 
been observed (23). However the reported p s (T,H) in YBCO by muon spin rotation appears 
not to exhibit the \[H behavior. A further study of p s (T, H) is clearly desirable. Again it 
may be more useful to introduce the scaling function by 



X(T ,H)- X ( T,0) 

l-21n2? + --- for ^ < l 



/(T/e) = AV - ; ' 7 = fi,(T/e) - 21n2(T/e) 



il n (1.622£)f forf»l (10) 
The scaling function I(T/e) versus T/e is shown in Fig.3. 

Nuclear spin lattice relaxation rate. - In the superclean limit T7 1 at low temperature is 
given by 



where 



2T J v iV ' V 2T' 

'TT A' 



(-x) 2j ( T A) (ii) 



J ( T / e ) = ^/ ^.9>)^( — s) (12) 
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Here T lr ^ is the nuclear spin lattice relaxation rate in the normal state. We have J(T/e) — > 1 

1 7T 2 T 

for T/e < 1 and J(T/e) -> ) 2 for T/e > 1. Therefore the nuclear spin lattice relaxation 

rate T-f^T, H)/T{^ takes - iJ behavior in the limit of T -> and - (J) 2 in the limit of 
T e Again it is more convenient to introduce the scaling function by 

- 1 ( r r tj\ _ r r- 1 ( r r r\\ 



IX 



G{T/e) ee ^ "^ [4 ;; 2 "'" ; = J(T/e) - -{T/ef 



~\ T^ + °((t) 2 ) forf »1 1 J 

We show the scaling function G(T/e) in Fig. 3. The detection of this scaling function will 
be very useful. 

Planar thermal conductivity. - Following fli| the thermal conductivity in the superclean 
limit (i.eT/A <C H/H C 2 < 1 ) within the conducting plane in the vortex state with the field 
configuration of H | c is obtained as 

K(T,H)/K n = — (Lf{^L) 2 + (lnf 4A ]) 2 } (14) 

where n n = is the thermal conductivity in the normal state, n the density of electrons 

m the electron mass, and T the scattering rate by impurity. In the limit of T — > 0, k(T, H)/n n 
Sv 2 

reduces to 2^ A 2 ^' ^ near nl ^ contrary to This linear field dependence as well as the 
quasi-scaling behavior of k/T 3 versus H/T 2 in Eq.(14) have not been seen in high-T c cuprates. 

On the other hand, a recent thermal conductivity data of Sr2Ru04 at T = 0.35 K exhibits 
a clear tf-linear dependence for H > 0.02 T Q. At a low magnetic field (H < 0.015 T) 
the thermal conductivity appears to be independent of H. The thermal conductivity in the 
superclean limit in the present model for T > e becomes almost independent of H and increases 
like /c~ T 3 ln(A/T). 

Further, this kind of behavior should be quite common to the unconventional supercon- 
ductors with the nodal superconductors (e.g. E 2u -state in UPt 3 [pi). In particular when 
e/T -C 1, k/T 3 decreases with increasing H almost linearly in H/T^. Indeed this behavior 
has been observed in the B-phase of UPt3 recently . 

The low temperature thermal conductivity in high-T c cuprates with the field H || c appears 
to be described by the one in the clean limit (i.e. H/ H C 2 <I/A< 1 ) rather than the one in 
the superclean limit. 

The thermal conductivity in the clean limit, on the other hand, is given by Jl4[ ; 

A 00 / v 2 eH , , /2A X , 4 A , * . „. 

^ = T( 1 + WA h ( 4 fc) h te» (15) 



Here Aoo is the order parameter in the absence of the impurity scattering at T = 0, koo = 
f Aoo"m ^ ne universal thermal conductivity in the limit of T — > and for H = 0. Indeed this 
~ H \n(const./ H) behavior is very consistent with the recent data from YBCO and Bi2212 
single crystals p6| ]. We show such a comparison in Fig. 4, where we took r/r c = 0.06, 0.188, 
and 0.33. Here T c = 0.8819 T c q is the critical scattering rate when the superconductivity 
disappears. 
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Fig. 4. - The thermal conductivity data in YBCO by Chiao et al. [26] are fitted with the the thermal 

conductivity Eq.(15) with T/r c = 0.06(— ), 0.188( ), 0.33(- • •), where □ for x = 0, O for x = 

0.006, A for x = 0.03 of YBa2(Cui_ x Zn ;!; )306.9. The H \n(const./ H) dependence describes very well 
the experimental data. 

Concluding Remarks. - Limiting ourselves in the configuration H || c (i.e. the magnetic 
field normal to the conducting plane) and in the superclean limit we obtain the expression of 
the thermodynamic quantities and Ty in NMR in vortex states in d-wave superconductors. 
Some of limiting behaviors for T — ► K have been well established in the vortex state of 
YBCO, though little work on the scaling behavior has been done. We have shown also the 
present model describes some features of the vortex state in Sr2Ru04, specially the presence 
of the nodal structure in the order parameter. This suggest that single crystals of Sr2Ru04 
will provide another testing ground of the present model. At the present, among three 2D 
/-wave states for Sr 2 Ru04, only one of them, A(k) ~ cos(c/c3)e ±l< ^, appears to be viable, since 
both the extremely small angle dependence of the upper critical field |3(| and the thermal 
conductivity |28, [H] with the field applied parallel to the conducting plane are incompatible to 
other two candidates. Also, we expect similar behaviors of vortex state in other unconventional 
superconductors [|l4| |2| . Therefore the exploration of the vortex state will bring further insight 
in unconventional superconductors. 
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